ID model for the dynamics and expansion of elongated Bose-Einstein condensates 



Pietro Massignanf] and Michele Modugncjj] 
INFM - LENS - Dipartimento di Fisica, Universita di Firenze 
Via Nello Carrara 1, 50019 Sesto Fiorentino, Italy 
(Dated: February 1, 2008) 

We present a ID effective model for the evolution of a cigar-shaped Bose-Einstein condensate in 
time dependent potentials whose radial component is harmonic. We apply this model to investigate 
the dynamics and expansion of condensates in ID optical lattices, by comparing our predictions 
with recent experimental data and theoretical results. We also discuss negative-mass effects which 
could be probed during the expansion of a condensate moving in an optical lattice. 
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I. INTRODUCTION 



The realization of Bose-Einstein condensates (BECs) 
in ID optical lattices provides an important tool to inves- 
tigate the superfluid and coherence properties of macro- 
scopic quantum systems. Recent experiments have al- 
lowed to probe static and dynamic features of these sys- 
tems such as the Josephson effect jl], ||, the emergence 
of squeezed states Q, the transition from superfluid to 
dissipative dynamics Bloch oscillations |J|, interfer- 
ence patterns produced by freely expanding condensates 
H fjj , and band spectroscopy |f| . 

These results have stimulated a great interest also from 
the theoretical point of view, and several models have 
been proposed to interpret the experimental data. In 
particular, new analytic results have been derived in the 
so called "tight binding regime" , based on the assumption 
that the condensate wave functions in different sites are 
well separated ||, §|. 

A very important point would be to compare these 
results with the numerical solution of the full 3D Gross- 
Pitaevskii equation (GPE) that describes the dynamics 
of BECs at very low temperatures, allowing also for the 
exploration of the weak binding regime. However, in gen- 
eral the solution of the 3D GPE is a formidable task that 
suffers severe limitations due to the large memory stor- 
age and computational times required to simulate the ex- 
perimental configurations. To overcome these difficulties 
the GPE is usually reduced to a ID equation by means 
of a static renormalization of the inter-atomic scattering 
length (l0[ [n], 12 , by using the fact that these systems 
are effectively one-dimensional. This allows for a realis- 
tic description of the trapped axial dynamics, but at the 
cost of loosing any information on the radial degrees of 
freedom. In general this reduction is indeed too rude, 
since the transverse dimensions actually play a crucial 
role both during the expansion, when the ID constraint is 
removed, and during the trapped dynamics, where the in- 
terplay between the axial and radial components is neces- 



sary to account for collective excitations, e.g. quadrupole 
oscillations. 

Recently the authors of |l3| have proposed an effec- 
tive ID wave-equation for BECs confined in cylindrically 
symmetric potentials, by using a different approach that 
keeps track also of the transverse dimensions of the con- 
densate. The derivation is based on a variational Ansatz 
and the wave function of the condensate is factorized in 
the product of an axial component and a gaussian ra- 
dial one. The axial wave function satisfies a ID non 
polynomial Schrodinger equation (NPSE) that gives re- 
markable results for the axial ground-state and dynamics 
of cigar-shaped condensates in cylindrically symmetric 
traps. However, the derivation in |l3| is based on the 
assumption of a static harmonic radial confinement, and 
therefore the NPSE is not suitable to study the dynam- 
ics under variation of the trapping potential, in particu- 
lar the ballistic expansion of the condensate after being 
released from the trap. 

In this paper we present a new model that combines 
the idea of gaussian factorization with the unitary scal- 
ing and gauge transformations of [|l4| allowing for 
a suitable description of both the transverse and longi- 
tudinal evolution of elongated BECs in the presence of 
a time-dependent harmonic potential plus an arbitrary 
axial component. The above transformations are used to 
reabsorb the evolution due to the time variations of the 
harmonic potential, the latter being replaced by an effec- 
tive confinement that makes feasible a gaussian factoriza- 
tion even in the case of a sudden release of the external 
confinement. The evolution of the system is described 
by an effective ID GPE that is dynamically rescaled (dr- 
GPE) by means of the transformations of |l4|, |lf|. De- 
spite the one dimensional character of this equation, the 
model is capable to account for the transverse dynamics 
through an algebraic equation for the radial width cou- 
pled to the rfr-GPE. Therefore this model represents an 
useful tool to describe both the dynamics and the expan- 
sion of BECs in elongated geometries. 

The Paper is organized as follows. In Section P we 
present the derivation of the rfr-GPE, and in Section [II A 
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we discuss its properties in the presence of a pure har- 
monic potential, by making a comparison with analyti- 
cal results in the Thomas- Fermi limit. In Section III B 
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we use the dr-GPE to discuss the dynamics and expan- 
sion of BECs in the presence of ID optical lattices, by 
comparing its predictions with recent experimental and 
theoretical results [^], ||, We also present a method 
to probe negative-mass effects during the expansion of a 
condensate in the presence of an optical lattice. 



II. THE MODEL 

We consider a condensate confined in a time-dependent 
harmonic trapping potential Uy lo with an additional axial 
component U\d 

U(r,t) = U ho (r,t) + U 1D (r z ,t) (1) 

= 2 mw i(*) r l + \ mu? ^) r l + Ui D (r z ,t). 

The Gross-Pitaevskii equation for the wave-function of 
the condensate can be derived from the functional [fill 



S[*] = 



dt <Fy 
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V 2 -U 



gN 
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where N is the number of condensed atoms, g = 
4irh 2 a/m the coupling strength, m the atomic mass and 
a the inter-atomic scattering length. 

In order to account for the time-dependent variations 
of the confining potential it is convenient to apply the 
unitary gauge and scaling transform of There- 
fore we introduce a rescaled wave function ^(x, t), de- 
pending on the rescaled spatial coordinates Xi = Tij \ (t) 
and related to the true wave function \I>(r, t) by 



*(r,t) 



.&Er^ *({ri/Aj(t)},t) 



(3) 



The scaling parameters Xj (t) are solutions of 



w?(0) 



A j (t)Ai(t)A 2 (t)A 3 (t) 



w?(t)Ai(t) (4) 



with the initial conditions A(0) = 1 and A(0) = 0. 

The use of this scaling equations has two advantages: 
(i) the evolution of \& due to the variation of the harmonic 
trapping potential is mostly absorbed by the scaling and 



gauge transform |lj (in the Thomas- Fermi limit |^| is 
frozen to its initial value) ; (ii) the rescaled wave function 
^> always evolves in the presence of a fictitious harmonic 
confinement, that depends on the trapping frequencies at 
t = 0, and therefore the gaussian factorization is applica- 
ble even in the case of ballistic expansion, when the real 
trapping is turned off. 

In case of elongated condensates we can obtain a ID 
effective model by factorizing the rescaled wave function 
\J/(cc, t) in the product of an axial component ip(z,£) and 
a gaussian radial component (f>(x, y,t;a(z, €)), allowing 
for an exact integration over the transverse coordinates 
©0 



(x,y;a(z,t)) 



with the normalization conditions 



dz\iP(z,t)f 



(5) 



(6) 



(7) 



With this choice the transverse size of \E' is characterized 
by the width a(z, t) that is expected to be a slowly vary- 
ing function of time, since most of the evolution due to 
the variation of the trapping potential is absorbed by the 
unitary transform (^) . We note that from Eq. (||) and (^) 
even the true wave function ^ can be expressed in the 
factorized form 

* (r,t) = <f>(r x ,r v ,t;Z(r z ,t))ip(r z ,t) (8) 

whose transverse width S is given by 

Z(r z ,t) = \ x (t)a(r z /\ z ,t). (9) 

The equations of motion for ip and er can be derived 
by using a variational procedure from the action (|^) , in- 
tegrating out the dependence on </>. The problem can 
be further simplified in the so called "slowly varying ap- 
proximation" , that is when the contribution given by the 
axial variations of the radial wave function <j> is negligi- 
ble, V 2 w Vj_(/> Jjl|. In this approximation the action 
becomes 



1 ^ 1 9N |^| 2 1 ^ 



2ma 2 X 2 j_ 2 A 2 A 2 L 2 2na 2 A Z A^ 



where it is possible to distinguish the contributions of the gaussian integration to the kinetic and potential energy 
(respectively the terms proportional to 1/er 2 and a 2 ). The requirement of stationarity of this functional yields the 
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equations for ip 
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+ U 1D (X z z,t) + 



\A 2 



mw^(O) 



raw, (0) 9 



gN 

27TCT 2 



^(*) (11) 



and for a 



10- 



a(z,t) = a ± A/A,(<) + 2a7V|V(^i)| 5 



(12) 



where a± is defined through the harmonic frequency at 
t = by a± = v / n/muj ± (0). From Eqs. @ and (||), the 
radial width E of the true wave function '5 is given by 



E(r„t) = a ± A ± (t)^A,(i) (1 + 2aiV|V"(r z , t)| 2 ). (13) 

By combining Eqs. ( pT| ) and ( |l2| ) we get a ID nonlinear 
Schrodinger equation for ip, depending on the parameters 
Xj solution of Eq. (^J) , that we call dynamically rescaled 
Gross-Pitaevskii equation (dr-GPE). This equation is en- 
ergy conserving and requires the same computational ef- 
fort of a simple 1D-GPE, since the numerical solution of 
Eq. (^) is straightforward. In the particular case of a 
time-independent harmonic potential Eqs. (|ll|) and Jl2] ) 
reduce to the NPSE of 0. 

The contribution due to the transverse part of the en- 
ergy, that has been neglected within the approximation 



V 2 
gral 



can be estimated by considering the inte- 



2m 



(|Vr 



2m A? 



2mA? 



dz 



(14) 

where er' = d z a. We have verified that, in the case of 
a general time-dependent harmonic trapping, the contri- 
bution of this term is a negligible fraction of the total en- 
ergy (less than 0.1%) in a wide range of trap anisotropics 
lu z /luj_ and number of atoms N. As we show in the next 
Section, even in the presence of an optical lattice along 
the axial direction, where the density can show deep mod- 
ulations on a short scale, the model gives quite accurate 
results. 

In the rest of the paper we consider some applications 
of the model and make comparisons with analytical re- 
sults in the Thomas-Fermi (TF) limit, full 3D numeri- 
cal solutions of the GPE and experimental results, where 
available. The parameters are chosen in the range of typi- 
cal experiments at LENS ^, 0, |^| . In particular, we model 
condensates formed by 2 • 1(T < N < 2 • 10 5 atoms of 87 Rb 
in a cigar-shaped configuration with axial and radial har- 
monic trapping frequencies respectively io z — 2ir ■ 9 Hz 
and oj± = 2ir ■ 92 Hz. 

To numerically solve the model, we use a Runge-Kutta 
algorithm for the scaling equations and a FFT split-step 
method for the dr-GPE evolution, mapping the wave 
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function on a discretized lattice 

state of the system is found by using a standard imagi- 
nary time evolution p6[ . 
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FIG. 1: Comparison of the dr-GPE and TF predictions for 
the radial size of a condensate in a pure harmonic potential 
as a function of the number of atoms N. The shift between 
the two curves is expected since the TF approximation trun- 
cates the tails of the wave function, leading to a systematic 
underestimate of the ground state size, especially for small N. 



III. APPLICATIONS 
A. Harmonic trapping 

We start by considering the case of a condensate con- 
fined in a pure harmonic potential. As discussed in p3[ , 
the ground state axial profile given by the gaussian fac- 
torization (^) reproduces the full 3D solution much bet- 
ter than other known ID model in both the weakly and 
strongly (TF) interacting regimes. The model also pro- 
vides reasonable results for the radial profile: the gaus- 
sian Ansatz, though slightly overestimating the central 
density, gives a prediction for the rms radius i? rTOS = 
\f(r\) really close to the analytical TF result even in 
case of very large number of atoms, as shown in Fig. |l|. 

A remarkable feature of the rfr-GPE is the possibil- 
ity to describe the interplay between longitudinal and 
transverse oscillations induced by modulations of both 
the axial and radial part of the trapping potential. To 
illustrate this aspect we consider the low-lying collective 
excitations induced by a sudden variation of the trans- 
verse external confinement. In particular we start with 
a condensate in the ground state of an anisotropic trap 
with uj±_ = 27r • 92 Hz and, at t = 0, we suddenly switch 
the transverse frequency to a final value of u±_ = 2ir ■ 80 
Hz, thus inducing shape oscillations involving both the 
radial and axial directions. For very elongated conden- 
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FIG. 2: Evolution of the axial (Z rm3 , top) and radial (R rm s, 
bottom) sizes of a condensate performing shape oscillations, 
respectively obtained solving the dr-GPE and the TF Eqs. (^) 
and in both graphs it is possible to observe the superpo- 
sition of the quadrupole and the faster transverse breathing 
frequencies (N = 2 ■ 10 s , u z = 2tt ■ 8.7 Hz). 



sates in the TF regime, these oscillations are character- 
ized by the quadrupole and transverse breathing frequen- 



cies, co 



Q - 



and 



,TB 



2lj±_ respectively ]16|. In 



Fig. |2| we show the behavior of th e rad ial and axial sizes 
(respectively R rms and Z rms = (z 2 ) ), comparing the 
dr-GPE predictions with the TF values. The agreement 
is remarkably good for both the frequency and amplitude 
of the oscillations; the slight shift between the two curves 
is expected, due to the truncation of the tails of the wave 
function in the TF approximation (see Fig.Jj]). Despite 
the ID (axial) character of the dr-GPE, Fig. pfshows that 
the equation well mimics the predicted coupling between 
the axial and radial degrees of freedom. We have checked 
that the same behavior can be reproduced by changing 
the axial trapping instead of the radial one, or by using 
a resonant drive. 

Another main difference between our model and any 
static ID renormalization of the GPE (1D-GPE) comes 
out when we consider the ballistic expansion of a con- 
densate suddenly released from the trapping potential. 
In fact, in the typical experimental regimes the free ex- 
pansion of a condensate is governed by the TF equations 

m 

^/<r?> T ,(t) = A i (t) A /<^7(0) (15) 

where Xi (i) are solutions of Eq. (Q) with ui z (t) — lj±_ (i) = 
0. In the case of a ID GPE a similar relation holds for 
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FIG. 3: Axial and radial sizes (respectively Z rms and Rrms) of 
the condensate during a free expansion from a pure harmonic 
potential, as obtained with the dr-GPE and within the TF 
approximation in ID and 3D. 



the axial size, except that now the only scaling parameter 
obeys 



X 1D {t) 



^ 2 (0) 



-col(t)\ XD {t). 



(16) 



One immediate consequence of Eq. (|T^) is that in the 
TF limit a ID GPE equation obtained through a sim- 
ple renormalization of the coupling constant, although 
suited to reproduce the (axial) ground state of the sys- 
tem, in general overestimates the repulsive character of 
the interaction during the expansion. This is shown in 
Fig. |^ where we plot the axial size of the expanding con- 
densate as obtained from the dr-GPE and the ID and 
3D TF scaling. The nice agreement between the dr-GPE 
solution and the full dimensional case is due to the use 
of the scaling equations and to the dynamically reduced 
nonlinearity of the equation ( p"l| ) that, in the strongly 
interacting (TF) limit aN\ip\ 2 3> X z , takes the form 



ihdtip(z) 



mco 2 z (0) 



A 2 A i 



z 2 + 



Ui D (X z z : t) 
3 



hco ± (0)V2aN\ii\ 



(17) 



The inset of Fig. || shows that even for the transverse 
dynamics, despite the ID character of the dr-GPE, the 
agreement with the 3D (TF) solution is remarkably good. 



B. Optical lattice 

We now turn to consider the case of a condensate in 
a ID optical lattice. Experimentally this is realized by 
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means of a retro-reflected laser beam which creates a pe- 
riodic potential of the form 

U 1D (r z ) = s-E r cos 2 (2^r z /A opt ) (18) 

where X op t is the wavelength of the laser, E r = 
h 2 /2mX 2 rpt is the recoil energy of an atom absorbing 
one lattice phonon and s is a dimensionless parameter 
controlling the intensity of the lattice. Here we chose 
A pt = 795 nm and intensities in the range < s < 6 as 
reported in the experiments j|, Q |) . 

The ground state of the system is characterized by 
deep periodic variations in the axial density and by pro- 
nounced modulations in the radial width, that are well 
reproduced by our model even for high values of lattice 
intensities. This is evident in Fig. || where we compare 
the axial density 

p(z) ee (\^) ± = I d 2 r ± |*| 2 (19) 

and the rms transverse radius, integrated only over the 
radial directions 



R± ms {z) ee y^iK =JJ d*r ± r±\*\* ; (20) 

with the full 3D solution. 

To address the issue of collective excitations we con- 
sider both dipole oscillations induced by a sudden dis- 
placement of the harmonic potential (here we use Az = 
5 ^m), and quadrupole modes excited as before by chang- 
ing the harmonic transverse confinement. 

Concerning the former, on general grounds one expects 
a frequency shift which can be explained in terms of a 
mass renormalization due to the s-dependence of the dis- 
persion relation in the lowest band of the periodic poten- 
tial j|, |2lJ : — > yjm/m*u)z '. In a recent paper it has 
been predicted, by using a tight binding Ansatz, that the 
same effective mass m* can account also for the modifi- 
cation of the quadrupole frequencies M: for elongated 
condensates in the TF regime one expects a linear re- 
lationship between the two frequencies, to® = yjb/2uj® , 
independently of the lattice intensity s. 

In Fig. U we compare the effective mass extracted from 
the dipole and quadrupole modes with that of a single 
particle in a periodic potential, near the Brillouin zone 
center ^J. This picture shows that the predicted rela- 
tion between and ui® is very well verified also in the 
weak binding regime of small s, and that the effects of 
mean field interaction and of harmonic trapping (both 
included in the dr-GPE description) give only a minor 
correction to the single particle picture in a uniform pe- 
riodic potential. Our results are also in good agreement 
with a recent experiment realized at LENS j^j , as shown 
in the inset of Fig. [|. 

Other interesting features of BECs confined in ID opti- 
cal lattices are the interference effects which take place af- 
ter the removal of the combined optical+harmonic trap- 
ping potential. In fact, the free expansion of the system 
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FIG. 4: Axial density p(z) (top) and rms transverse radius, 
averaged only on the transverse coordinates and plotted as 
function of the axial coordinate (bottom); the insets show the 
central region of the condensate (N = 5 ■ 10 4 , s — 5). 
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FIG. 5: Effective mass m* , normalized to the unperturbed 
value m, as a function of the lattice intensity s, obtained re- 
spectively from the dipole (continuous line) and quadrupole 
(dashed) frequencies, and compared to the case of a single 
particle in a periodic potential (dotted). In the inset the fre- 
quency of the dipole mode is compared with the experimental 
data and theoretical curve (dashed-dotted) from |H|. 
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FIG. 6: Axial density of the central peak, as obtained respec- 
tively with the dr-GPE, a Schrodinger propagation and the 
1D-GPE introduced in pj| , compared with the experimental 
distribution of || (the wings are due to the presence of a small 
thermal component), after t exp — 29.5 ms. 




FIG. 7: Structure of the lowest bands for a single particle in 
a ID periodic potential with s = 2 (solid lines), compared to 
the free particle case (dashed). The inset shows the behavior 
of the effective mass m*, inversely proportional to the band 
curvature, in the first band for s — 2. 



causes the single condensates forming the initial array 
to overlap and leads to the formation of neat interfer- 
ence patterns, showing up as lateral peaks moving with 
quantized speeds v n = ±n • 2h/X op tm (with n integer) 
||. This is essentially a ID effect, which reflects the mo- 
mentum distribution of the ground state along the axial 
direction, and it is well reproduced by our model, re- 
garding both the position and the population of the first 
lateral peaks. An interesting point concerns the axial 
expansion of the central peak. A careful analysis of the 
expanded density profile shows that, compared to other 
existing ID equations, the dr-GPE gives a better esti- 
mate of the axial size of the central peak. In fact, as 
can be seen in Fig. ^, a (free) Schrodinger propagation of 
the ground state underestimates the axial width, while 
the statically renormalized 1D-GPE introduced in p| , 
though appropriate to describe the ground state, largely 
overestimates the axial size after the expansion, as dis- 
cussed in the previous section. This is indeed due to 
the fact that the meanfield term of the dr-GPE has a 
linear character (cx |V>|, see Eq. [It]) which lies between 
the Schrodinger (cx and the GPE (cx \tp\ 2 ) cases. 

Our model gives accurate predictions also for the radial 
expansion of the central peak, which is essentially gov- 
erned by the TF scaling solution, in agreement with the 
experimental values and the theory of ||. 

We conclude this section by discussing a time-of-flight 
method that could be used to probe the dispersion rela- 
tion in the presence of a periodic potential P3[ , allowing 
for the investigation of negative-mass effects during the 
expansion in an optical lattice. The idea is the follow- 
ing: we start with a condensate in the ground state of 
a pure harmonic potential, and we accelerate it up to a 
given velocity v — Tiq/m in the axial direction, e.g. by 
exciting a dipole motion in the harmonic potential. Then 
we switch off the trap and adiabatically ramp an optical 



lattice along the longitudinal direction up to a desired 
height s. Here we consider the case of a linear ramp of 
15 ms, which is sufficient to ensure the adiabaticity of the 
process pi| for the velocities considered here. We have 
indeed verified that applying the reverse ramp after the 
first one, the additional Fourier components of the wave 
function at q ± 2lqs (with I integer) disappear. 

With this procedure we can load the condensate in a 
Bloch state of axial quasimomentum hq and band index 
n, and therefore explore effects of effective mass (which 
depends on s, n and q) far from the Brillouin zone center, 
during the expansion of the system. In particular, since 
meanfield effects are ruled out after few milliseconds, the 
radial and axial expansion are almost decoupled, and the 
latter can be accounted for by the single particle mass 
renormalization, as shown in Fig. ^ || |TJ. The advan- 
tage of loading the condensate in the lattice only after 
the acceleration procedure is that in the reverse case the 
occurrence of dynamical instability effects prevents the 
possibility of accelerating the condensate above a critical 
velocity v c < vb 22|, where vb — TiqB/m = h/mX is the 



Bragg velocity. 

In Fig. [s] we show the evolution of the axial size of the 
condensate during the expansion in the optical lattice, for 
various initial velocities v. The behavior of the system 
is essentially determined by the single particle effective 
mass which characterizes the diffusive (kinetic) term in 
the evolution equation (as stated before, meanfield ef- 
fects become negligible after few milliseconds of expan- 
sion). Indeed, the figure shows that when q approaches 
the value q — 0.6895 qs at which the single particle ef- 
fective mass for s = 2 becomes infinite (see Fig. 0), the 
expansion in the direction of the lattice gets frozen. For 
values between q and qs, m* becomes negative, and the 
condensate starts contracting. On the contrary, when 
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FIG. 8: Evolution of the axial width of a condensate which is 
loaded during the expansion in a ID optical lattice, for various 
velocities: q/qs = 1.3, 1.2, (continuous line), 0.4, 0.685 (dot- 
ted), 0.75,0.8, respectively from top to bottom. The lattice 
intensity is ramped linearly in the first 15 ms of expansion, 
and then kept fixed at the value s = 2 for other 10 ms (see 
the inset). 
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FIG. 9: Aspect ratio after 25 ms of expansion, for a con- 
densate loaded adiabatically in a ID optical lattice of height 
s = 2, as a function of the condensate velocity. 



pect ratio after a fixed evolution time. Indeed, as shown 
in Fig. [I the aspect ratio is characterized by a finite 
jump across the first Brillouin zone boundary, which re- 
flects the discontinuity in the band curvature (the ef- 
fective mass passes from a negative to a positive value). 
With the same technique one could probe the band struc- 
ture also for higher values of q and/or s. In the latter 
case a longer ramp time could be required in order to 
ensure an adiabatic loading |p3. 



IV. CONCLUSIONS 

We have presented an effective ID model for the evolu- 
tion of an elongated Bosc-Einstein condensate in time de- 
pendent potentials whose radial component is harmonic, 
in the presence of an arbitrary axial potential. The 
model exploits the scaling and gauge transformations of 
|l4], [l^] , that reabsorb most of the evolution due to the 
time-dependent variation of the harmonic potential. This 
makes possible a gaussian factorization for the radial part 
of the condensate wave function, even in case of a free ex- 
pansion after the release from the trap. 

The evolution of the system is described by a non- 
linear ID wave equation (<ir-GPE) for the axial wave 
function, dynamically rescaled by means of the aforemen- 
tioned transformations. Despite the ID character of the 
rfr-GPE, the model also accounts for the radial dynam- 
ics of the system, that is necessary to correctly describe 
collective oscillations and the free expansion of the con- 
densate. 

The accuracy of the model has been tested by compar- 
ing its predictions with the TF solution in case of pure 
harmonic trapping and with recent experimental and the- 
oretical results for the dynamics and the expansion of 
BECs in ID optical lattices. We have also shown that 
the quasimomentum dependence of the effective mass has 
dramatic consequences, which could be easily accessed 
in the experiments by accelerating and letting expand a 
condensate in an optical lattice. 
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the condensate is loaded in the second band, that is for 
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